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Abstract 


'^Thls  report  deals  with  the  use  of  many  sensors  grouped  into  two  or  more 
clusters  (subarrays)  for  the  determination  of  bearing  and  range  to  an  acoustic 
source  producing  sinusoidal  or  narrowband  Gaussian  signals.  The  spacing 
between  subarrays  is  typically  very  large  compared  with  the  acoustic  wavelength, 
but  the  spacing  between  sensors  in  any  one  subarray  is  comparable  with  or 
smaller  than  the  acoustic  wavelength.  The  noise  is  Gaussian  and  spatially 
incoherent .■^^^Che” following  are  some  of  the  more  general  conclusions. 

A.  Sinusoidal  signals 

1.  When  only  two  subarrays  are  used  the  estimator  which  optimally  combines 
the  outputs  of  all  sensors  offers  little  advantage  over  a system  which 
obtains  separate  bearing  estimates  from  each  subarray  and  determines 
source  location  by  simple  trlangulatlon. 

2.  When  three  or  more  subarrays  are  available  an  estimator  which  forms 
conventional  beams  at  each  subarray  and  combines  these  coherently  comes 
close  to  the  optimum.  Hence  there  is  little  or  no  incentive  to  build 

a much  more  complicated  system  using  individual  sensor  data  as  Inputs. 

Each  subarray  acts,  in  effect,  as  a single  sensor  with  liq>roved  signal 
to  noise  ratio.  The  directional  properties  of  these  equivalent  sensors 
are  of  only  minor  Importance. 

3.  The  accuracy  of  a system  employing  three  or  more  subarrays  could,  in 
principle,  be  far  greater  than  that  of  a single  subarray.  In  practice,  the 
factor  discuased  in  the  next  paragraph  severly  limits  this  potential 
Improvement. 

4.  Whan  the  aignal  is  a pure  sinusoid  estimation  accuracy  is  limited 

by  aablgulty  problems.  These  are  caused  by  the  fact  that 


differential  distances  to  a pair  of  sensors  are  Indistinguishable  If 
they  differ  by  a multiple  of  the  signal  wavelength.  In  the  absence  of 
source  motion  the  only  good  way  to  resolve  such  ambiguities  appears  to 
be  the  use  of  signals  which  are  not  single  frequency  sinusoids. 

B.  Narrowband  Gaussian  Signals 

1.  When  the  spacing  between  subarrays  Is  large  compared  with  the  average 
wavelength  of  the  signal,  most  of  the  useful  Information  appears  to  be 
contained  In  the  envelope  of  the  received  signal.  Only  the  envelopes 

of  the  subarray  outputs  were  therefore  used  In  coherent  estimation 
procedures . 

2.  Formal  evaluation  of  the  optimal  estimator  working  with  envelope  Infor- 
mation Is  a very  difficult  problem.  However,  lower  and  upper  bounds  on 
Its  mean  square  error  can  be  calculated,  the  former  by  working  with  the 
pre-envelopes,  the  latter  by  evaluating  a particular  Instrumentation. 

3.  One  of  the  key  questions  Is  whether  the  coherent  combination  of 
envelopes  from  several  subarrays  yields  a better  estimate  than  conventional 
procedures  using  the  sensors  In  only  subarray,  but  now  utilizing  carrier 

as  well  as  sideband  Information.  The  answer  Is  straightforward  for  range 
estimation:  The  combination  of  envelopes  from  several  subarrays  yields 
a far  more  accurate  estimate  under  most  reasonable  conditions.  The  case 
of  bearing  estimation  Is  less  clear  cut.  It  Is  possible  to  find  entirely 
reasonable  combinations  of  parameters  for  which  either  one  or  the  other 
scheaw  will  be  preferable.  For  the  single  array  the  far  field  bearing  accu- 
racy Is  detcntlned  primarily  by  the  ratio  of  array  length  to  average 
signal  wavelength.  The  corresponding  quantity  for  the  system  eaq>loying 
severe!  suberreys  Is  the  ratio  of  spacing  between  subarrays  to  wavelength 


of  the  highest  modulation  frequency 


4.  K version  of  the  ordinary  (2  element)  split  beam  tracker  was  used 

to  set  upper  bounds  on  the  minimum  mean  square  error  In  delay  estimation. 
This  tracker  employs  the  squared  envelopes  of  two  subarray  outputs  as 
Its  Input  signals.  For  very  large  tlme-bandwldth  products  and  large 
Input  signal  to  noise  ratios;  the  upper  and  lower  bounds  differ  by  about 
6 db.  Since  the  lower  bound  Is  almost  certainly  not  realizable,  this 
suggests  that  the  split  beam  envelope  tracker  Is  not  too  far  from  the 
optimum  In  this  Important  range  of  parameters.  When  the  signal  to 
noise  ratio  at  each  subarray  (after  beamforming)  falls  below  unity  the 
tracker  performance  deteriorates  rapidly  compared  with  the  lower  bound. 

5.  When  the  tlme-bandwldth  product  Is  only  moderately  large,  the  split 
beam  tracker  error  approaches  a finite  minimum  as  the  signal  to  noise 
ratio  Increases.  It  comes  close  to  this  minimum  when  the  signal  to  noise 
ratio  exceeds  the  tlme-bandwldth  product.  The  lower  bound  predicts  that 
the  mean  square  error  approaches  zero  as  the  signal  to  noise  ratio  grows 
Indefinitely.  In  this  case,  therefore,  the  split  beam  tracker  suffers 

a definite  performance  loss  relative  to  the  lower  bound  when  the  signal 
to  noise  ratio  becomes  large. 

6.  An  Improved  (balanced)  version  of  the  split  beam  tracker  overcomes 
the  problem  mentioned  In  5.  For  signal  to  noise  ratios  well  In  excess 

of  0 db  Its  mean  square  error  differs  from  the  lower  bound  by  an  essentially 

/ 

constant  factor  of  about  6 db. 


1.  Introduction 

This  report  summarizes  work  carried  out  under  contract  N66001-75-C-0210 
between  Yale  University  and  the  Naval  Undersea  Center,  San  Diego.  Major 
portions  of  the  material  are  adapted  from  a progress  report  and  three  letter 
reports  submitted  earlier. 

The  study  Is  concerned  with  the  localization  of  an  acoustic  signal  source 
by  means  of  an  array  consisting  of  many  sensors  grouped  Into  two  or  more 
subarrays.  The  acoustic  source  radiates  a signal  which  consists  either  of 
pure  alnuAolds  or  of  very  narrowband  noise- like  waveshapes.  The  wavelength 
of  the  signal  sinusoid  (average  wavelength  In  the  case  of  a narrowband 
signal)  Is  likely  to  be  small  compared  with  the  dimensions  of  any  one  subarray 
but  Is  probably  comparable  to  (or  larger  than)  the  spacing  between  Individual 
sensors  In  the  subarray.  Thus  the  subarrays  have  beam  patterns  with  signi- 
ficant directional  properties.  The  spacing  between  subarrays  Is  large 
coiq)ared  with  the  dimensions  of  any  one  subarray  and  therefore  extremely 
large  coiq>ared  with  the  acoustic  wavelength.  The  frequency  of  the  signal 
sinusoid  (or  center-frequency  of  the  narrowband  signal)  Is  assumed  to  be 
known  and  fixed  during  the  observation  Interval.  Thus  there  Is  no  effort 
to  determine,  or  exploit,  any  effects  of  source  motion.  Information  concerning 
source  location  (bearing  and  range)  Is  contained  only  In  the  relative  delay 
of  the  signal  observed  at  various  sensors. 

All  of  the  analysis  Is  carried  out  under  the  assumption  that  the  signal 
wavefront  la  perfectly  coherent  from  subarray  to  subarray  (and  hence  certainly 
over  any  one  aubarray) . To  a first  approximation  one  could  think  of  a less 
than  perfectly  coherent  signal  as  being  composed  of  two  parts:  One  which 
la  perfectly  coherent  and  one  which  la  completely  Incoherent.  The  Incoherent 
part  simply  adds  to  the  noise,  so  that  the  effect  of  partial  Incoherence  Is 


n 

• 
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merely  a reduction  in  signal  to  noise  ratio. 

Only  the  most  primitive  noise  model  Is  used:  The  noise  at  each  sensor 
Is  Gaussian,  and  the  noise  conqponents  received  by  different  sensors  arc  statis- 
tically Independent.  For  the  analysis  of  sinusoidal  signals  the  noise  at 
each  sensor  Is  further  assumed  to  be  white.  This  assumption  Is  probably 
quite  realistic:  Practical  sea  noise  spectra  should  be  very  nearly  flat  over 
the  narrow  frequency  bands  of  Interest.  Spatial  Incoherence  Is  a more  doiibt- 
ful  assumption,  at  least  within  any  one  subarray.  However,  the  advantages  to 
be  gained  by  exploiting  such  short  range  spatial  structure  of  the  noise  are  both 
limited  and  well  understood.  To  Include  this  effect  In  the  present  analysis 
would  greatly  complicate  the  results  without  shedding  much  additional  light 
on  the  questions  of  primary  Interest. 

The  basic  objective  of  the  study  Is  to  determine  what  lnq>rovements  In  bearing 
and  range  accuracy  can  be  obtained  by  combining  the  information  supplied  by  various 
subarrays.  The  required  combination  can  be  accomplished  at  several  different 
levels: 

(1)  The  most  general  approach  Is  to  regard  all  sensor  outputs  as  potential 
Inputs  Into  a single  optimal  estimator.  If  the  number  of  sensors  is  at 
all  large,  the  resulting  Instrumentation  becomes  extremely  complex. 

(2)  At  the  other  extreme,  one  can  use  conventional  beamforming  at  each 
subarray  to  obtain  separate  bearing  estimates.  These  can  then  be  combined 
by  purely  geometrical  techniques  to  obtain  an  estimate  of  source  location. 

(3)  As  an  Intermediate  procedure,  one  can  beamform  at  each  subarray,  and 
then  use  coherent  combinations  of  the  beemformer  outputs  to  obtain  the 
estimate  of  source  location.  This  approach  effectively  treats  each  sub- 
array as  a single  (directional)  sensor  In  a larger  array  whose  output 

is  now  processed  In  optimal  fashion.  If  the  nuaber  of  subarrays  Is  not 


very  large,  the  resulting  instrumentation  might  not  be  excessively 
complex. 
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All  of  these  procedures  have  been  studied  to  a greater  or  lesser  extent. 

We  shall  begin  this  presentation  with  the  most  general  approach  (1).  It 
will  quickly  become  apparent,  however,  that  Its  theoretical  advantages 
over  (3)  are  quite  marginal  and  certainly  Insufficient  to  Justify  the  vastly 
greater  complexity.  The  real  choice  Is  therefore  between  (2)  and  (3).  The 
decision  between  these  two  Is  much  less  clearcut.  The  major  advantage  of 

the  coherent  combination  of  array  outputs  envisioned  in  (3)  is  that  it  exploits  j 

I 

the  large  baseline  provided  by  the  subarray  spacing.  On  the  other  hand,  the 

fact  that  the  subarray  spacing  is  likely  to  be  very  large  compared  with  the  { 

acoustic  wavelength  leads  to  ambiguities  (multiple  lobes  In  the  effective  i 

beam  pattern)  which  are  very  difficult  to  resolve.  In  the  absence  of  source 

motion,  signal  bandwidth  appears  to  furnish  the  only  visible  mechanism  for 

dealing  with  this  problem.  The  tradeoff  between  (2)  and  (3)  can  therefore 

be  stated  very  roughly  as  follows: 

(a)  One  can  work  with  arrays  of  relatively  small  dimensions  (the  subarrays)  ^ 

using  a wavelength  determined  by  the  center  frequency  of  the  signal 

''i'; 

spectrun. 

(b)  One  can  exploit  the  full  spacing  of  the  subarrays,  but  the  relevant 
wavelength  Is  now  determined  by  the  highest  modulation  frequency  and  Is 
therefore  much  larger  than  the  wavelength  associated  with  the  center 
frequency. 

For  different  sets  of  parameter  values  the  advantage  can  shift  from  one  to  the 
other  of  theae  procedures.  A more  detailed  examination  of  these  questions  Is 
a primary  ala  of  this  report. 
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2.  Basic  Theory  and  Deflaitlons 


source 


Xh 


o- 


subarray  2. 


su  bar 


Figure  1 

The  geometrical  relations  between  source  and  receiver  locations  are 
Indicated  In  Fig.  1.  We  are  Interested  In  estimating  the  source  bearing 
e and  range  r relative  to  some  conveniently  chosen  origin  of  coordinates 
0.  The  available  data  are  the  outputs  of  M sensors,  x^(t) . . .x^(t) . Each 
of  these  consists  of  a signal  component  and  a noise  component  and  Is 
observed  In  the  time  Interval  (0,  T) . In  line  with  the  earlier  comments 
the  signal  components  at  all  sensors  are  asstawd  to  be  Identical  In  form 
except  for  delays  determined  by  strictly  geometrical  considerations.  We 
do  allow  different  signal  attenatlons  to  the  various  subarrays  and 
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can  therefore  discuss  the  effect  of  using  subarrays  with  different  signal 
to  noise  ratios. 

There  Is  always  a degree  of  arbitrariness  In  the  definition  of  an 
"optimum"  estimator.  We  shall  require  at  least  that  the  estimator  Intro- 
duce no  systematic  error,  so  that  the  actual  error  will  tend  to  zero  as 
the  observation  time  T Increases  without  bound.  In. statistical  terms, 
we  require  that  our  estimator  be  "unbiased".  Beyond  that,  we  take  the 
usual  approach  of  characterizing  estimator  performance  by  a mean  square 
error  criterion.  This  enables  us  to  use  the  well-known  Cram%i>Rao 
technique  for  obtaining  the  desired  bounds  on  estimation  error. 

Specifically,  let  each  sensor  output  be  represented  by  a set  of  n 
numbers.  These  might  be  the  n time  8aiiq>les  x^(t^) . . .x^(t^) , a set  of  n 
Fourier  coefficients,  or  any  other  convenient  representation.  Let  ic  be 
the  complete  data  vector  of  dimension  M«n,  consisting  of  n saiq>les  from 
each  of  the  M sensors.  The  likelihood  function  A Is  defined. by 

A - p(x| e,  r)  (1) 

Thus  A Is  the  probability  density  of  the  received  data  vector  x when  the 
true  values  of  the  unknown  parameters  are  6 and  r.  We  must  now  evaluate 
the  Fisher  Information  matrix 


£{}  stands  for  the  expectation  of  the  bracketed  quantity. 

The  Cramir-Rao  theory  asserts  that  the  mean  square  estimation  errors 

2 * 2 “ 

D (6)  of  6 and  D (r)  of  r have  absolute  lower  bounds  given  by 
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D^(0)  > (3) 

> [3‘^]22 

The  symbol  [ Indicates  the  element  In  the  1th  row  and  jth  column 
of  the  matrix  In  brackets. 

The  lower  bounds  are  entirely  general.  Independent  of  the  particular 
estimation  procedure  used,  assuming  only  that  the  estimate  Is  unbiased. 

The  element  [J  Inverse  Fisher  matrix  measures  the  correlation 

between  the  bearing  and  range  estimates. 

One  cannot  assert  In  general  that  the  lower  bounds  of  Eqs.  (3)  and  (4) 
are  necessarily  attainable.  The  general  theory  Indicates  that  for  large 
observation  times  T,  the  (realizable)  maximum  likelihood  estimator  achieves 
mean  square  errors  approaching  those  of  the  lower  bounds.  "Large  T"  means. 

In  practice,  a time  large  compared  with  the  correlation  time  of  the  noise  (and 
the  signal,  If  It  Is  random).  This  condition  Is  quite  likely  to  be  satisfied 
In  situations  of  practical  Interest.  In  this  sense  the  Cram^r-Rao  bound  can 
therefore  be  regarded  as  a meaningful  limit  which  can  at  least  be  approached 
with  Instrumentations  of  sufficient  complexity.  There  Is,  however,  one 
Important  qualification:  Most  practical  instrumentations,  even  when  they 
purport  to  Implement  the  maximum  likelihood  estimator,  do  not  locate  the 
global  maximum  of  the  likelihood  function.  They  find  a local  maximum,  a 
point  at  which  the  derivative  of  the  likelihood  function  with  respect  to 
e and  r la  zero.  When  the  signal  Is  very  narrowband,  there  are  many  local 
maxima.  In  using  Cramer-Rao  as  a realizable  bound  we  are  therefore  assuming 
that  we  have  successfully  resolved  the  ambiguity  problem  posed  by  multiple 
peaks  of  comparable  height.  As  the  signal  bandwidth  diminishes,  the  difference 


Id  height  of .Adjacent  peaks  decreases  and  the  resolution  problem  becomes 
more  and  more  difficult. 


Sinusoidal  Signal.  Two  Subarrays 


As  a first  specific  case  we  consider  a source  radiating  a purely 
sinusoidal  signal  ^ of  known  frequency  but  unknown  phase 

s(t)  - Aq  slnCUpt  - ij>).  (5 

The  receiver  consists  of  two  subarrays  which  will  be  taken  as  line  arrays 
with  equally  spaced  elements  In  order  to  simplify  trigonometric  manipulations 
Source  bearing  and  range  are  measured  relative  to  a system  of  coordinates 
centered  at  the  left  end  of  the  first  subarray,  as  Indicated  In  Fig.  2. 


>c:  source 


Figure  2. 


In  the  light  of  the  discussion  Just  completed  this  problem  has  a rather  arti- 
ficial flavor.  With  a purely  sinusoidal  signal  there  Is  no  practical  way  to 
resolve  the  ambiguities.  To  calculate  how  accurate  a system  would  be  which  did 
resolve  them  may  appear  to  be  an  academic  exercise.  The  reason  for  pursuing  the 
argument  nevertheless  Is  that  It  exhibits  In  a particularly  simple  setting  a 
number  of  problems  which  remain,  but  become  much  less  transparent,  when  the 
signal  Is  sufficiently  complex  to  permit  Identification  of  the  global  peek  In 
the  likelihood  function.  These  problems  Include  the  coupling  between  bearing 
and  range  eatlmatlon  and  the  merits  of  coherently  combining  beam  outputs  as 
opposed  to  individual  sensor  outputs. 
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We  assume  Chat  the  source  range  from  each  subarray  (r,  p)  Is  large 
compared  with  the  subarray  dimensions  (L^,  ^2^*  ^ ^ result  the  angles 

6 and  do  not  vary  significantly  over  any  one  subarray.  An  equivalent 
observation  is  that  the  signal  wavefront  is  essentially  planar  over  any 
one  subarray.  We  do  not  assume  that  the  spacing  D between  subarrays 
is  necessarily  small  compared  with  r or  p. 

If  we  designate  Che  signal  travel  time  from  the  source  to  the  1^^ 
sensor  by  the  waveshape  received  at  that  sensor  will  be 

x^(t)  - slnlu>Q(t  - T^)  - (>1  + n^(t)  (6) 

n^(t)  is  Che  Gaussian  white  noise  associated  with  the  1^^  sensor  and  is, 
by  assumption.  Independent  of  all  other  noise  components.  Only  the  delay 
depends  on  r and  6. 

The  phase  angle  ^ of  the  transmitted  signal  is  without  intrinsic 
Interest,  but  it  is  essential  to  the  computation  for  Che  following  reason: 

If  ^ is  taken  as  zero  (or  as  some  other  known  quantity) , the  time  of  zero 
crossings  at  the  source  is  known  and  one  can  determine  the  range  to  any  one 
receiver  by  observing  the  time  of  zero  crossings  at  that  point  (to  multiples 
of  a wavelength) . Ranges  to  two  receivers  would  locate  the  source  in  bearing 
and  range.  In  practice,  this  is  clearly  nonsense.  One  does  not  know  the 
time  of  zero  crossings  st  the  source  and  observations  at  two  sensors  can 
only  locoes  the  source  on  the  hyperbola  corresponding  to  the  observed 
dlffersncial  delay. 

There  are  two  basic  approaches  to  the  problem  Just  outlined.  One  can 
regard  ^ aa  another  fixed  but  unknown  parameter  to  be  estlaiated  in  addition 
to  S and  r.  This  enlarges  the  Fisher  information  matrix  [Bq.  (2)]  to  3 x 3 
format  and  conaidarably  complicates  the  Inversion  process.  Alternatively 
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one  can  treat  ^ aa  a random  parameter  [perhaps  uniformly  distributed  over 
(0,  2ir)]  and  average  over  its  distribution.  This  preserves  the  2x2 
dimension  of  the  information  matrix,  but  leads  to  exceedingly  cumbersome 
Integrals  in  its  individual  elements.  The  first  approach  is  chosen  as 
the  lesser  evil. 

To  Illustrate  the  computation  of  the  elements  of  J,  consider  the  term 


gj3_log_Aj.  take  the  noise  spectral  level  at  the  i^**  sensor  as  N and 
assume  a noise  white  over  a frequency  band  W.  Then  Myquist  rate  samples 
taken  at  Intervals  of  1/2  W are  statistically  independent  and,  because  of 
the  assumed  spatial  Incoherence  of  the  noise,  the  likelihood  function 
assumes  the  simple  form 

M n 1 2 

A - C exp  {-  E E (t.)  - s (t  )]  } (7) 

1-1  j-i  1 J ^ J 

where  C is  a normalizing  constant.  s^(t)  is  the  signal  received  at  the 
1^^  sensor  and  is  given  by 


*l(t)  - \ sin[u>Q(t  - T^)  - ♦].  (8) 

Differentiating  the  natural  logarithm  of  (7)  twice  with  respect  to  0 
one  obtains 


38.  (t.) 


3V(t,) 


3 log  A _ 1*  r _L.  / r_-i_JL_i2  a.  r \ ill 

30^  1-1  J-1  “r  1 J 1 J 3q2 

Since  the  noise  has  zero  mean  the  last  term  in  Eq.  (9)  disappears  when 

the  expected  value  is  taken  and  one  is  left  with 


(9) 


(10) 


Substituting  Eq.  (8)  into  Eq.  (10)  one  obtains 
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2 2 

.2,  . M n A.  0).  , 3t.  « 

- 1 - I E cos2(«  (t  - , ) - ♦)  (U) 

36^  i-1  j-1  i w j i 


For  large  W there  are  many  time  samples  In  each  signal  period  and  one  can 
1 


approximate  the  J sum  by  an  Integral. 

2 2 

«2,._  . n M A_,  u- 


2 2 T 

M A u.  3t  2 r 2 

“ ' 1 ^ 

^ 0 

In  practice  the  observation  time  T would  almost  certainly  cover 
many  periods  of  the  signal  sinusoid.  Then  one  can  approximate  to  a high 
degree  of  accuracy^ 
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The  other  elements  of  the  Fisher  Information  matrix  are  computed  In 
entirely  analogous  fashion.  The  result  Is 
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Note  that  the  time  Increment  from  sample  to  sample  Is  dt  - 1/(2W). 
The  result  is  exact  when  T is  a multiple  of  the  signal  period. 


A tedious  but  basically  straigJit f orward  computation  now  yields  the 
determinant  of  J and  the  various  cofactors  required  for  the  evaluation 
of  the  Inverse  matrix.  The  terms  of  primary  Interest  are 
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(15) 

(16) 

(17) 

(18) 


stands  for  the  cofactor  of  the  matrix  element  in  the  rth  row  and 
•Cth  column. 

Eqs.  (15)  - (18)  are  quite  general.  They  do  not  use  the  particular 
array  configuration  of  Fig.  2,  nor  do  they  require  the  source  to  be  remote 
compared  with  the  subarray  dimensions.  We  now  introduce  these  assumptions 
in  order  to  put  the  various  partial  derivatives  into  a relatively  simple 
form. 

When  the  index  J is  associated  with  the  first  subarray  in  Fig.  2 we 
have  the  approximation 


T 


i 


Jd  cos  6) 


(19) 


where  c is  the  velocity  of  sound  and  d the  distance  between  adjacent 
elements.  The  required  partial  derivatives  are  wherefore 
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(20) 


J I sin  e 


(21) 


Vfhen  the  Index  k is  associated  with  the  second  subarray  the  same 
approximation  yields 


- -(p  - kd  cos  \|i) 


(22) 


p and  ^ are  source  range  and  bearing  from  the  second  subarray  as  Indicated 
in  Fig.  2.  The  required  partlals  are  therefore 


_ ^^k  3p  ^"^k 

3r  **  3p  3r  3^  3 


r ° if) 


(23) 


. !!k  iP  + . i/3£.  + kd  sin  >P  ^) 

36  3p  36  ^ 3i;i  36  c^as  36^  * 


(24) 


Simple  trigonometric  manipulations  based  on  Fig.  2 yield  explicit 
expressions  for  the  various  partial  derivatives  ^ in  Eqs.  (23)  and  (24) 
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cos  a 


- cos  a 


- r sin  o 
36 


(25) 

(26) 

(27) 

(28) 


a Is  the  angular  separation  of  the  subarrays  as  seen  from  the  source  [see 
Fig.  2]. 

The  next  step  is  to  substitute  Eqs.  (23)  - (28)  into  Eqs.  (15)  - (18) 

^ They  are,  of  course,  simply  Che  elemenCs  of  the  Jscoblsn  describing  the 
cransformsClon  from  the  (r,  6)  coordlnsce  system  to  the  (p,  \|i)  coordinate 

system. 


««  Jti. ' if ' 


I-,- 


and  use  the  results  to  evaluate  Eqs.  (3)  and  (4).  Considerable  simpli- 
fication can  be  achieved  by  assuming  that  the  signal  to  noise  ratio  is 
constant  over  any  one  subarray.  Specifically  we  shall  write 

A 2 r 

i _ j -jjp"  when  1 is  in  subarray  1 

Ni  ; 1 

A,^ 

V.  —jJ”  when  i is  in  subarray  2 

It  is  also  helpful  to  assume  that  the  numbers  of  elements  in  the  two 
subarrays,  and  M2  respectively,  are  much  larger  than  unity.  Then 
one  can  approximate 


Mid  - L^, 


where  and  L2  are  the  array  lengths. 

Even  with  these  approximations  the  required  algebraic  manipulations 
are  still  extremely  tedious  and  only  the  final  results  are  therefore 


reproduced  here. 
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The  newly  introduced  symbol  A stands  for  the  wavelength  of  the  signal 
sinusoid. 

Insplte  of  their  coinplexity  Eqs.  (31)  and  (32)  yield  some  immediate 
checks  against  physical  Intuition: 

1)  If  ■*  0,  l.e.  the  second  subarray  Is  not  receiving  any  signal, 

Eq.  (31)  reduces  to  the  very  simple  expression 

2 “ 12 

D^(e)  > ^ r (33) 

2 ^ ^1  2 

2 ^ 

The  combination  ^ measure  of  the  post-beamformlng  signal 

to  noise  ratio.  With  that' Interpretation  Eq.  (33)  has  all'  of  the 
features  one  would  expect  to  find  In  the  minimum  mean  square  bearing 
error  attainable  with  a linear  array. 

2)  When  A^  0 the  right  side  of  Eq.  (31)  tends  to  Infinity.  This 
conclusion  becomes  reasonable  as  soon  as  one  recalls  that  6 was  defined 
as  the  source  bearing  relative  to  subarray  1.  The  second  siibarray  pro- 

' vldes  a measurement  of  ^ which  may  be  quite  accurate.  It  does  not  per- 
mit a range  measurement  and  hence  does  almost  nothing  to  locate  the 
source  In  bearing  from  the  position  of  subarray  1. 

3)  The  right  side  of  Eq.  (32)  tends  to  Infinity  when  either  A^  or 
A2  tend  to  zero.  Since  the  wavefront  curvature  over  each  subarray 
Is  nagllgible,  one  subarray  alone  cannot  produce  a range  estimate. 

In  practlca  both  subarrays  would  have  to  receive  substantial  signal 
levels  if  the  overall  systesi  is  to  perform  slgnlflcsntly  better  than  one 
of  its  components.  For  simplicity  we  consider  the  special  case 
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r- 


(34) 


S/N  therefore  stands  for  the  post-beamformlng  signal  to  noise  ratio  of  each 

$ 

subarray  and  W Is  the  processed  frequency  band.  If  we  Ignore  terms  of  the 
order  L^/p  and  1>^/^  we  can  simplify  Eqs.  (31)  and  (32)  drastically: 


0^(0)  > 
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(35) 


2 <?  2 "2  r 2 2 
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D^(r)  ^ 


1 - cos  a 


2 2 
I ^ T ^ 2 

2 S 2 ^2  r 7 2 

4ir  TW  ^ sin  e + 2 2 ^ J 
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1 + cos  g 
1 - cos  a 


2 2 * 

rV(0) 


(36) 


The  first  Interesting  comparison  Is  between  Eq.  (35)  and  Eq.  (33). 

If  and  r = p the  addition  of  the  second  subarray,  even  when  used 

In  optimal  fashion,  produces  only  a modest  reduction  In  bearing  error. 

One  gains  the  obvious  advantages  from  having  one  subarray  or  the  other 
oriented  more  nearly  broadside  with  respect  to  the  source,  but  there  Is 
no  apparent  benefit  from  the  enormous  Increase  In  overall  array  dimension 
to  the  order  of  D rather  than  L^. 

Eq.  (36)  tells  an  analogous  story.  Given  two  similar  subarrays  at 
similar  distances  from  the  source  one  could  simply  use  each  one  separately 
to  obtain  a bearing  estimate  and  then  locate  the  source  as  indicated  in 
Figure  3.  Treating  the  sub array outputs  coherently  appears  to  yield  only 
modest  benefits  relative  to  this  simple  geoisetrlcal  coad)ination  of  bearing 
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estimates . 


Figure  3. 

It  Is  not  difficult  to  account  for  the  failure  of  the  large  overall 
array  dimension  to  Influence  the  range  and  bearing  accuracy  more  drastically. 
One  way  to  approach  the  problem  Is  to  consider  the  correlation  between  the 
bearing  and  range  estimates.  The  matrix  J ^ la  the  covariance  matrix  of 
the  best  Joint  estimate  of  0,  r and  It  la  therefore  a straightforward, 
though  once  again  tedious,  matter  to  compute  the  correlation  coefficient 

A A 

between  the  aatlmatlon  errora  of  6 and  r. 


P 


Or 


'12 


'^*^11'^22 


(37) 


Computations  have  been  carried  out  only  for  aqual  algnal  to  nolac  ratios 
at  the  two  arrays,  lu  %iblch  case  one  finds  a correlation  coefficient 
differing  from  unity  only  by  tema  of  the  order  very 
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small  quantities.  One  can  Interpret  these  observations  as  follows: 

Each  of  the  subarrays,  even  In  the  optimal  Instrumentation,  produces 
an  effect  not  very  different  from  that  of  a single  directional  sensor. 

By  combining  the  array  outputs  coherently  one  establishes,  to  a high 
degree  of  precision,  a hyperbolic  locus  of  constant  differential  distance 
from  the  two  subarrays  (solid  line  In  Fig.  4).  The  segment  of  the  hyperbola 


Figure  4 

on  which  the  source  la  actually  located  Is  Identified  only  by  the  directional 
properties  of  each  "sensor",  i.e.  by  the  estimation  accuracy  of  the  separate 
subarrays  (dotted  lines  In  Fig.  4).  With  coherent  processing,  therefore, 
the  area  of  uncertainty  Is  a very  narrow  strip  along  the  hyperbola  and 
within  the  dotted  quadrilateral.^  With  separate  processing  of  the  subarray 

outputs  It  Is  the  entire  dotted  quadrilateral.  The  hyperbola  traverses. 

^ We  are,  once  again.  Ignoring  the  ambiguity  problem,  which  will  be  examined 
In  section  S. 


within  Che  quadrilateral,  a variety  of  bearing  and  range  values  not 
dramatically  smaller  than  Chose  covered  by  the  entire  quadrilateral. 
This  accounts  for  the  limited  Improvement  In  bearing  and  range  accuracy 
due  to  coherent  processing. 
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4.  Sinusoidal  Signal.  Three  Subarrays. 

If  two  subarrays  Identify  a hyperbola  of  probable  source  location, 
three  subarrays  should  define  two  such  hyperbolas  and  their  intersection 
should  pinpoint  the  source  location  to  a high  degree  of  precision.  It 
appears  promising,  therefore,  to  examine  the  configuration  shown  in  Fig.  5. 
As  in  Eq.  (35)  we  imnediately  discard  terms  of  order  L^^/r  and  L^/p.  We 
also  ignore  terms  of  the  order  L^/X  and  L^/X  as  negligible  compared  with 
Dj^/X  and 


Figure  5. 

With  equal  post-beamformlng  signal  to  noise  ratio  S/N  at  each  subarray  we 
than  obtain 


2 * 

D^(0)  > 


12[(coa  g - 1)^  cos(B  - 1)^  (cos  a - cos  g)^l 


(38) 


2 S 2 


^ sin  Yj^(co8  o - 1)  - sin  Y2(co8  6-1) 


Comparison  with  Eq.  (35)  shows  the  expected  result:  The  accuracy  is  now 


22 


governed  by  the  ratio  (D^  -t*  rather  than  the  much  smaller  quantities 

p*rallel  observations  ean  be'made'ooncemlag  the 

range  accuracy. 

Eq.  (38)  goes  a long  way  towards  resolving  a question  raised  much  earlier 
In  this  report:  How  should  one  combine  the  Information  provided  by  the  v.irious 
subarraya  In  order  to  achieve  a reasonable  compromise  between  accuracy  and 
complexity  [see  discussion  on  p.  2]?  In  Eq.  (38)  the  subarrays  contribute 
to  bearing  accuracy  only  by  enhancing  the  poat-beam  forming  signal  to  noise 
ratio.  One  would  obtain  exactly  the  same  result  by  considering  an  array 
composed  of  only  three  omnidirectional  sensors,  each  with  signal  to  noise 
ratio  S/N.  The  absolute  optimum  estimator,  which  combines  all  sensor  outputs 
In  optimal  fashion,  has  a mean  square  error  differing  from  the  above  only 
by  terms  which  we  discarded  as  negligible  In  arriving  at  Eq.  (38).  It  is 
clear,  therefore,  that  very  little  Is  gained  by  using  the  subarrays  Co  do 
anything  mors  sophisticated  than  to  enhance  signal  to  noise  ratio  by 
conventional  beamforming. ^ 


It  appears  plausible  to  extrapolate  this  statement  to  noise  fields  which 
are  not  apatially  incoherent  over  one  subarray.  One  should  now  construct  a 
rsaaonabla  approximation  to  tha  optimum  beamformsr  for  the  subarray  and  then 
eoiAina  subarray  contributions  as  before. 


On  the  exle  Joining  the  two  subarrays  the  distance  between  adjacent 
hyperbolas  Is  X/2.  It  remains  of  that  order  until  r Is  much  larger  than 
D,  then  It  begins  to  grnw  linearly  with  r.  The  number  of  possible  hyperbolas 
within  the  bearing  defined  uncertainty  region  (shaded  area  of  Fig.  3)  would 
presumably  be  large  and  would  not  diminish  with  Increasing  distance  from 
the  subarrays.  Source  location  Is  therefore  not  constrained  to  the  neighbor- 
hood of  a single  line,  as  suggested  by  Fig.  4,  but  to  neighborhoods  of  many 
near-parallel  lines.  The  utility  of  the  Information  gained  by  coherent 
processing  of  the  two  array  outputs  therefore  becomes  very  marginal. 

The  situation  Is  not  Improved  materially  through  the  addition  of  a 
third  subarray.  We  now  have  two  Intersecting  families  of  hyperbolas,  yielding 
a large  number  of  points  as  possible  source  locations.  Ambiguity  therefore 
appears  to  be  one  of  the  real  limiting  factors,  at  least  In  the  problem 
setting  considered  thus  far. 

The  obvious  means  of  avoiding,  or  at  least  reducing,  ambiguity  Is  to 
work  with  signals  containing  two  or  more  frequency  components  not  harmoni- 
cally related.  If  the  signal  consists  of  several  pure  sinusoids,  the  analy- 
sis of  sections  3 and  4 can  be  used  with  relatively  minor  extensions. 
s^(t)  Is  now  a sum  of  sinusoids  which  must  be  substituted  into  Eq.  (10). 

Since  each  of  these  sinusoids  has  the  same  delay  the  only  change  In 

2 2 

Eq.  (13)  will  be  replacement  of  T by  a more  complicated  term  In- 

volving amplitudes  and  frequencies  of  all  signal  components.  The  remain- 
der of  the  computations  would  not  change  materially. 

If  the  signal  does  not  consist  of  fixed  sinusoids,  but  covers  a con- 
tinuous band  of  frequencies,  the  ambiguities  can  be  resolved  formally 
In  the  limit  of  extremely  long  observation  times.  In  order  to  gain  some 
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Insight  Into  the  seriousness  of  the  problem  for  realistic  observation  times 
we  consider  the  following  very  simple  case. 

We  model  the  signal  as  a narrowband  Gaussian  process  of  center  frequency 
u)^  and  bandwidth  Au.  The  receivers  are  two  subarrays,  modelled  as  omnidirec- 
tional sensors  with  the  signal  to  noise  ratio  achieved  after  conventional 
beamforming  We  are  Interested  In  determining  the  differentiated  delay 
of  the  signal  between  the  two  sensors. 

For  observation  times  large  compared  with  the  correlation  times  of  signal 

2 

and  noise  this  estimation  problem  Is  solved  optimally  by  the  ordinary  two 
element  split  beam  tracker  sketched  In  Fig.  7.  H Is  a filter  matched  to 
the  spectral  properties  of  signal  and  noise  (Eckart  filter). 


di  Cato  r 


Figure  7 


^ A more  realistic  model  would  work  with  directional  sensors,  but  this  should 
have  little  influence  on  the  eehiguity  problem  discussed  here. 

^ W,  J,  Bangs,  Array  Processing  with  Generalised  Beamformers,  Report  Ho.  1, 
Contract  M0001A0-72-C-1293,  NUSC  Sept.  1971. 


26 


When  the  delay  Ax  Introduced  at  the  output  of  sensor  2 matches  the 
actual  delay  Ax  of  the  signal  at  sensor  1 relative  to  that  at  sensor  2, 

A 

the  expected  output  E(z)  Is  zero.  Thus  one  need  only  adjust  Ax  until  the 
output  reaches  a null  and  use  Its  value  as  an  estimate  of  Ax  and  hence  of 
differential  delay. 

If  the  signal  Is  a pure  sinusoid^ E(z) , considered  as  a function  of 

A 

Ax , will  simply  be  a sinusoid.  It  will  then  have  an  Infinite  number  of 

A 

nulls  at  all  values  of  Ax  corresponding  to  differential  delays  differing 
from  the  true  value  by  a multiple  of  the  signal  period.  If  the  signal 
Is  a narrowband  process  one  would  expect  basically  similar  behavior,  with 

A 

one  Important  difference:  The  plot  of  E(z)  vs  Ax  will  now  have  an 
envelope  as  suggested  In  Fig.  8.  It  will  still  have  a series  of  nulls. 


Figure  8 

but  the  fact  that  the  peak  closest  to  the  correct  null  Is  higher  than  all 
others  provides  a means  of  resolving  the  ambiguity. 

If  the  signal  bandwidth  Is  very  small  compared  to  Its  center  frequency, 
the  envelope  of  Fig.  8 decays  slowly  and  the  difference  In  height  of 
successive  maxima  becomes  difficult  to  detect.  As  a measure  of  this 
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W.  J.  Bangs,  p.  101,  p.  157. 

Tha  filter  U in  Fig.  7 la  here  simply  a perfect  bandpass  filter  matched  to 
the  signal  band. 

3 * 

The  exact  expression  for  arbitrary  dr  -dt  Is  not  difficult  to  obtain,  but  Is 

quite  cumbersome. 


difficulty  we  compute  the  difference  in  the  values  of  E(z)  at  two  successive 
maxima  and  divide  It  by  the  "typical  fluctuation  noise",  the  rms  value  of  the 
fluctuation  of  z. 

An  analytical  expression  for  E(z)  can  be  extracted  from  the  previously 

cited  reference^.  Assuming  a signal  spectrum  of  level  S flat  over  the  band 

2 

Wq  - du  ^ u ^ Uq  + do)  we  obtain 


“2  sin  du)(dT  -At)  cos  <*>q(At  -At) 


[Aw  cos  <*>q(At  -dT)cos  Aw(At  - At) 


~ “QSin  “q(At  -At)  sin  Aw(dT  -At)] 


For  small  Aw/w^  the  first  two  maxima  of  E(z)  (above  At  - At)  should  occur 
very  close  to  tt/(2Wq)  and  5it/(2Wq)  respectively.  Evaluating  E(z)  at 
these  two  points  and  using  the  approximation  Aw/w^  « 1 one  obtains 


^ 2 ,Aw,3 

”"o  • 


'At-At  - 


At-At 


The  standard  deviation  of  z should  not  vary  greatly  over  an  interval 

of  one,  or  even  several,  cycles  in  Pig.  8.  We  can  therefore  work  with  its 
* 3 

value  at  At  •At.  The  required  analytical  formulation  Is  again  contained 
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r 


A 

ft 


In  the  cited  report  . To  the  same  order  of  approximation  as  In  Eq.  (40) 

»/  * 2 f «o 

where  N Is  the  white  noise  power  level  at  each  receiver. 

Hence  the  figure  of  merit  for  our  ability  to  resolve  the  ambiguity 
Is 


E(z) 


- E(z) 


At  - At 


2m, 


At  - At« 


5-n 

2ii>_ 


D(z) 


m,. 


(42) 


If  we  are  to  distinguish  the  height  of  successive  peaks  and  hence  Identify 
the  proper  null,  we  must  certainly  have  a figure  of  merit  well  In  excess 
of  unity.  For  a 12  bandwidth  and  high  signal  to  noise  ratio  the  required 
condition  Is 


TAw  > 


10^  a 3.2  X 10^ 
3 S S/N 

N 


(43) 


Unless  Uq  (and  hence  Am)  Is  large  this  can  easily  lead  to  exceedingly  large 
2 

observation  times. 

It  Is,  of  course,  not  necessary  to  resolve  adl scent  nulls  of  Fig.  8 
^ W.  J.  Bangs,  p.  38. 

2 

The  strictly  linear  S/N  dependence  for  large  signal  to  noise  ratios  results 
from  the  use  of  D(z)  at  At  > At.  A more  exact  coaq>utation  would  result  In  a 
figure  of  merit  which  saturates  at  vary  high  S/N. 
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in  order  to  gain  some  benefit  from  the  coherent  processing  procedure. 
When  the  signal  frequency  Is  high  and  differential  distances  of  a wave- 
length are  therefore  quite  small,  considerable  gains  can  be  made  by 
Identifying  the  location  of  the  true  null  within  10  or  20  periods  of 
Fig.  8.  The  extent  to  which  this  Is  possible  clearly  depends  on  the 
envelope  of  Fig.  8.  What  we  have  considered  here  Is  only  one  possi- 
ble way  of  utilizing  the  envelope  information. 
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6.  Narrowband  Signals;  General  Considerations* 

For  the  remainder  of  this  report  we  shall  assume  that  the  signal  Is 
a narrowband  Gaussian  process  whose  spectrum  Is  centered  about  some 
frequency  large  compared  with  the  bandwidth.  If  the  bandwidth  Is 
extremely  small  (smaller  than  the  Inverse  of  the  observation  time)  the 
signal  looks  essentially  like  a fixed  amplitude  sinusoid  throughout  the 
observation  Interval.  The  problem  Is  then  essentially  the  same  as  that 
treated  In  the  previous  sections,  the  only  difference  being  that  the 
amplitude  of  the  received  sinusoid  Is  Itself  a random  variable. 

If  a sinusoid  of  known  amplitude  Is  unsatisfactory  (because  of  ambiguity) , 
a sinusoid  of  the  same  frequency  but  unknown  amplitude  will  certainly 
be  no  Improvement.  It  follows  that  the  only  bandwldths  of  practical 
interest  are  those  which  are  appreciably  larger  than  the  Inverse  of  the 
observation  time. 

On  the  other  hand,  It  Is  necessary  to  keep  In  mind  that  signals  of 
practical  Interest  may  have  a bandwidth  which  Is  quite  small  compared  with 
the  center  frequency  ui^.  Unless  is  very  large,  or  the  observation  time 
Is  unrealistically  long,  the  time -bandwidth  products  will  therefore  not 
exceed  unity  by  many  orders  of  magnitude.  We  shall  find  that  this  leads  to 
the  breakdown  of  an  analytical  assumption  which  has  been  made  very  widely 
In  the  study  of  optimal  detectors  and  estimators. 

We  now  turn  to  the  role  of  the  signal  envelope.  We  have  seen  that  the 
narrowband  estimation  problem  consists  of  two  distinct  components 

a)  Separation  of  the  global  maximum  of  the  likelihood  function  from  other 

suxima. 

b)  Estimation  of  the  coordinates  of  the  global  maxlmua. 
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The  argument  of  Section  5 suggests  that  errors  due  to  a)  will  far  outweigh 
those  due  to  b)  for  typical  parameter  values.  With  a few  widely  separated 
sensors  (subarrays)  the  ability  to  separate  maxima  of  the  likelihood  function  Is 
governed  by  the  envelope  of  the  signal  correlation.  [See  Fig.  (8).]  This 
envelope.  In  turn,  la  generated  by  the  signal  envelope.  One  Is  therefore 
led  to  the  conjecture  that  one  can  envelope  detect  the  beamformed  outputs 
of  each  subarray  and  generate  the  location  estimate  from  a coherent  combina- 
tion of  the  envelopes  without  sacrificing  any  Information  of  real  value. 

An  alternative  statement  of  this  point  of  view  Is  the  following:  It  Is 
hopeless  to  atteiq>t  measurement  of  differential  distances  to  one  wavelength, 
or  even  a few  wavelengths.  The  only  really  relevant  Information  In  the  signal 
Is  therefore  contained  In  Its  slowly  varying  envelope.  We  shall  adopt  this 
point  of  view  and  assume  from  here  on  that  each  subarray  output  has  been 
envelope  detected  before  further  processing.  We  have  no  precise  measure  of 
what  has  been  thrown  away  In  doing  this,  but  we  suspect  strongly  that  the 
sacrifice  Is  small. 

Under  these  ground  rules , the  most  general  estimator  assumes  the  form 
shown  In  Fig.  9.  The  designer  Is  given  the  envelopes  as  Inputs  and 
must  use  them  to  obtain  the  best  estimates  of  bearing  and  range.  Since  the 
signal  Is  known  to  have  narrow  bandwidth.  It  Is  reasonable  to  assume  that  prior 
to  envelope  detection  each  beamformer  output  has  been  processed  by  a bandpass 
filter  cantered  about  and  of  a bandwidth  at  least  comparable  with  that 

of  tha  signal.  Thus  the  noise  as  well  as  the  signal  at  the  envelope  detector 
Inputs  has  narrowband  properties. 


i 

* 

i 


33 


if  Bignal  and  iioiue  uro  iiidopciidoiiL  (lauHaian  procuHaea  the  boamformur 
outputs  x^(t)  are  Gaussian  processes  which  can  be  written  In  envelope  and 
angle  form  as  follows^ 


X.  (t)  - X (t)  cos  u-t  + X (t)  sin  u)„t 
1 c^  0 . 0 


where 


and 


- R^(t)  cosLoj^t  + 6^(t)] 


- /x  ^(t)  + 


(t)  + X ‘(t) 
^1  ®1 


*8 

e^(t)  . arctan 


(44) 


(45) 


(46) 


X (t)  and  X (t)  are  Gaussian  random  processes,  each  with  a spectrum  obtained 
1 ®1 

from  Che  spectrum  of  Che  original  narrowband  process  x^(t)  by  downshifting  It 


(Jq  rad/sec. 


It  Is  apparent  from  Eq.  (45)  that  R^(c)  Is  not  a Gaussian  random  process. 

In  fact,  It  Is  a simple  matter  to  demonstrate  that  for  any  given  t the  random 
variable  R^(t)  has  a Rayleigh  distribution.  To  characterize  the  estimator  Input 
statistically  one  needs  the  Joint  distributions  of  Che  R^(tj)  for  many  obser- 
vation Instants  t^.  It  Is  not  difficult  to  write  this  distribution  down 
formally,  but  It  Is  extremely  difficult  Co  Infer  from  It  Che  performance  of 
the  maxlimtm  likelihood  estimator  or  even  Che  Carmdr-Rso  bound.  Efforts  Co 
proceed  along  these  lines  were  entirely  unsuccessful.  An  alternative  approach 
was  therafore  taken. 

If  one  assumes  that  the  pre-envelope  processes  x (t)  and  z (t),  1 ••  1,  2,  3, 


Sea,  for  Instance,  S.  0.  Rice,  MaChemaclcal  Analysis  of  Random  Noise, 

B.  S.  T.  J.,  Jan.  1975. 
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are  separately  available  at  the  estimator  Input  (which  they  are  not) , the  data 

vector  becomes  Gaussian  and  the  estimation  procedure  Is  well  understood.  One 

can  Immediately  write  down  the  relevant  estimation  errors  from  results  available 

In  the  literature.  Since  the  Input  processes  R^(t)  can  be  constructed  from 

the  X (t)  and  x (t)  [by  using  Eq.  (45)],  but  not  vice  versa,  such  calculations 
‘^l  ®1 

yield  lower  hounds  on  the  errors  attainable  with  the  Instrumentation  of  Fig.  9. 

Upper  bounds  on  estimation  error  can  be  obtained  by  considering  any  parti- 
cular# physically  realizable  Instrumentation.  If,  with  careful  choice  and 
some  luck,  one  arrives  at  an  Instrumentation  whose  errors  are  not  too  far  above 
the  lower  bounds,  then  one  can  claim  to  have  bracketed  the  estimation  errors 
attainable  with  Fig.  9 quite  accurately.  At  the  same  time  one  has  described 
an  Instrumentation  which  comes  resbonably  close  to  the  optimum.  That  this 
Is  Indeed  possible,  at  least  over  Important  ranges  of  the  Interesting  para- 
meters, Is  demonstrated  in  a later  section. 


7«  Lower  Bound  on  Estimation  Error 


Bangs  has  solved  the  problem  of  optimum  bearing  and  range  estimation 

using  a set  of  sensors  each  of  which  puts  out  a Gaussian  random  process^. 

To  adapt  this  result  for  our  purposed  we  need  only  recognize  that  each 

sensor  In  our  lower  bound  approximation  produces  two  Gaussian  random 

processes,  x (t)  and  x (t).  We  must  therefore  determine  the  cross-correlation 
*^1  ®1 

between  these  two  processes.  A straightforward  con4>utatlon  yields 

00 

K{*c  (t^)}  ■ I G^(w)  8ln[(u)  - “Q)(t2  “ t^)]  dm  . (47) 

^ ^ 0 

G^(w)  Is  the  spectrum  of  the  narrowband  process  x^(t)  produced  by  the  1^^ 

sensor.  When  t^  " t^^  the  Integrand  of  Eq.  (47)  Is  zero,  confirming  the 

well-known  results  that  the  two  pre-envelope  functions  x (t)  and  x (t) 

*^1  ®1 

are  statistically  Independent  when  observed  at  the  same  time.  A similar 
claim  Is  not  true  in  general  for  t^  5*  t^.  However,  If  the  spectrum  of  the 
narrowband  process  Is  symmetrical  about  l.e. 


G^(w)  - G(u)  - Wq) 


(48) 


with  G(x)  an  even  function,  then 


E{x  (t.)x  (t  ) } 

Cl  1 s^  2 


G(x)  8ln[(t2  - t^)  x]  dx  = 0 


(49) 


The  only  approximation  Implied  In  Eq.  (49)  la  that  the  bandwidth  of  the 

process  Is  sufficiently  small  compared  with  so  that  the  lower  limit 

can  be  replaced  with  In  that  case  the  x (t)  yield  a data  set  completely 

®1 

Independent  of  the  x (t).  One  can  then  work  with  either  member  of  the 

*^1 

pair  and  reduce  the  calculated  mean  square  errors  by  a factor  of  two.  The 
2 

results  of  Bangi  aaelysls  are  now  applicable  with  only  trivial  modification. 

^ U.  J.  Bangs,  Array  Processing  with  Generalized  Beamformers,  Report  No.  1, 
Contract  N00140-72-C-1293,  NUSC,  Sept.  1971. 

^ W.  J.  Bangs,  p.  97. 
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a)  Two  subarrays 


In  Che  simplest  case  there  are  only  two  subarrays.  Since  we  are  modelling 

2 

each  of  these  as  an  omnidirectional  sensor  we  can  only  estimate  differential 
delay  (l.e.  we  can  only  Identify  the  hyperbola  on  which  the  source  Is  located). 
The  general  expression  for  Che  estimation  error  Is 


r,\h  ■ I 

0 


-1 

du)  } 


(50) 


S(w)  and  N(u)  are  the  signal  and  noise  spectra  of  the  pre-envelopes. 

Their  shape  Is  that  of  the  narrowband  spectra  downshifted  In  frequency  by  (Jq. 

If  the  signal  to  noise  ratio  Is  constant  (S/N)  over  Che  signal  band 
0 _<  (0  Aw  and  zero  elsewhere  a simple  calculation  yields 


2 “ 
D^(t) 


3'it 

2 


1_ 

(Aw) 


1_ 

2 TAw 


1 -h  2S/N 
(S/N)^ 


(51) 


In  order  to  gain  some  Intuitive  feeling  for  Che  magnitude  of  this  error, 
suppose  Chat  the  source  Is  remote  compared  with  Che  subarray  spacing,  so 
Chat  error  In  differential  delay  can  be  translated  Into  bearing  error. 


2 '' 

D^(e) 


-J.  _i_  ( 1 2S/N 

Sir  TAw  '■D  sin  0^  (S/N)^ 


(52) 


Here  Is  Che  wavelength  of  the  highest  modulation  frequency,  related  to 
Aw  by 


X 


M 


2nc 

Aw 


(53) 


4 

As  an  example,  suppose  D « 10  A,  X ■ 1000  X (O.IZ  bandwidth)  and  S/N  - 10. 

n 

with  these  numerical  values  the  rms  bearing  error  Is 


0 0 

0(0) ■ degrees  (54) 

•'^Aw  sin  0 

Thus,  even  with  this  extremely  narrowband  signal,  one  does  not  require 

excessively  large  TAw  products  In  order  to  obtain  reasonable  bearing  accuracy. 

2 ’ ———————————————————— 

With  a signal  to  nolae  ratio  equal  Co  the  posC-beamformlng  signal  to  noise 
ratio  of  Che  aubarray. 
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The  matrices  and  are  composed  of  purely  geometrical  factors. 


8Ti  3t 

^^^IJ  “ 30  " 30 


9Tj^  3t 


^^Ij  “ 3r  ” 3r 


Is  the  travel  time  of  the  signal  from  the  source  to  the  1 subarray. 

\ 

The  parameter  In  Eqs.  (55)  and  (56)  describes  the  dependence  between 
the  bearing  and  range  measurements.  It  Is  given  by  the  expression 

[TrCl^^)]^ 

Pre  " t r ^59) 

Tr(T^^)Tr(T^^ 

■ 0 Indicates  Independence  of  the  measurements, 

P^g  ■ 1 Implies  complete  dependence. 

The  matrix  elements  given  by  Eqs.  (57)  and  (58)  are  easily  computed. 

They  assume  a very  simple  form  when  expressed  In  terms  of  the  angles  a and  3 
subtended  at  the  source  by  the  various  array  pairs. 


2 r . 

•5^  - - am  a. 


^^3  r 


^^1  1 ^^2  1 

“57  “ c*  “j;  " 7 “• 


®^3  1 . 

-5^  - - cos  3. 


Straightforward  computations  now  yield  the  geometrical  factors  required 
In  Eqs.  (55)  and  (56) 


Tr(TQTQ^)  ■ {sin^o  + sin^g  + (sin  a - sin  3)^)  (62) 

c 

Tr(T^.T^^)  • — j {(1  - cos  a)^  + (1  - cos  3)^  + (coso  - cos  3)^}  (63) 

c 
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p _ {sin  a(co8  g - 1)  + sin  g(cos  g - 1)  -<•  (sin  g - sin  a)  (cos  8 - cos  a)] 
[sin^a  + sin^P  + (slna  -sin  3)^][(1-cob  a)^  + (l-cos  6)^  + (cosa  -cos  3)^] 


(64) 


A numerical  example  Is  worked  out  In  Fig.  11  for  S/N  - 10.  Note  that  the 
angles  a and  3 decrease  from  the  top  to  the  bottom  of  the  table. 


! a 

3 

^r0 

~ D(0) 

O 

60 

120“ 

0.57 

0.13 

3.1“ 

45“ 

90“ 

0.80 

0.28 

16“ 

e 

O 

45“ 

0.92 

1.6 

36“ 

10“ 

20“ 

0.92 

7.2 

75“ 

5“ 

10“ 

0.92 

29 

151“ 

Fig.  11 

This  implies  that  the  source  is  more  and  more  distant  compared  to  the  spacing 

between  subarrays.  The  correlation  coefficient  Increases  at  first, 

but  then  quickly  approaches  an  asymptotic  value  as  the  angles  a and  3 

decrease.  For  a remote  source  the  correlation  coefficient  is  quite  close  f 

to  unity,  indicating  strong  coupling  between  bearing  and  range  measurements.^ 

The  fourth  column  of  Fig.  11  gives  the  normalized  range  error  (in  units 
of  modulation  wavelength) . Once  the  source  is  remote  compared  with  the 
Interarray  spacing  (last  two  entries)  the  range  error  increases  with  the  ji 

square  of  the  actual  range  as  one  would  expect.  The  last  column  gives  the  B 


It  is  possible  to  reduce  the  coupling  by  choosing  the  reference  point  (from 
which  bearing  and  range  are  measured)  near  the  center  of  the  overall  array 
rather  than  at  one  end. 


* 


f 


normalized  bearing  error.  Note  chat  Che  normalization  Includes  Che  value 
of  r.  For  large  ranges  Che  actual  bearing  error  becomes  Independent 
of  r,  as  one  would  anticipate. 

The  analytical  expressions  of  Eqs.  (55)  and  (56)  are  quite  complex. 
Very  much  simpler  forms  result  If  the  actual  range  happens  to  be  large 
compared  with  Che  Interarray  spacing  (although  no  assumption  to  that  effect 
Is  made  In  obtaining  the  estimate).  ^ 


2 * 
D''(0) 


=.  3 1 

16it  TAu) 


1 3 S/N  /m  . 2 1 ^ 

(S/N)^  ‘*13  ^ ■ ^r© 


sln^(0  -!(')  + sln^O  - sin  0 sln(0  - ^) 

di3  d^3 


(65) 


2 * 

D‘'(r)  =■ 


_3  JL_ 
2ir  TAu 


1 3 S/N  v2 

(S/N)-^  d-^^^ 


i 


(66) 


A **12  A A 2 **12  2 2 2 

sin**  (0  - Ip)  + (■T^)  sin  0 + [sin  (0  - ip)  - (“J^)  sin  0] 

**13  **13 


Suppose,  In  particular,  that  the  3 subarrays  are  arranged  In  a straight  line 
(ip  ■ 0)  and  are  equally  spaced  (d^^  ■ 2 d^^^)  • ^^en 


d2(0)  . L±Am  (_:h ^^2  _JL 

(S/N)^  **13  ® ^ ^ 


r0 


.2,*.  12  1 It  3 S/K  , Si'  .2  2 1 

° <'>  ■ te ~2  < 2 , 2 ) ' Trr:  • 

(S/M)  d^3  sin  0 r0 

For  large  actual  ranges  the  correlation  coefficient  p . reduces  to 

ro 


(67) 

(68) 
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{8in^(0-i/»)+(-T^)^8ln^0  - ■^(8in(0-^)  + 8lnQ]  8ln(0-i(»)  8in0}^ 
^13 I ^13  *^13 

[8in^f0-ip)+(-j^)^8ln^0  - 8ln(0-i<i)8in6  ] [8in^(0-ip)+(-T^)^8in^0-  (-r^)  ^8ln^(0-^)8in^0] 

o ®i  o d-  o Q«  ^ 


This  l8  Independent  of  r and  only  a function  of  array  geometry.  With  ^ > 0 and 
d^2  * *2  *^23  aS8ume8  the  numerical  value 

p - — . (70' 

It  may  be  Intereatlng  to  compare  the  bearing  eatlmate  obtained  by  the  3 
8ubarray  configuration  with  that  product  by  the  2 subarray  configuration  (under 
a far- field  aasumptlon)  . We  therefore  conalder  the  aame  combination  of  para- 

4 

meter  valuea  as  In  the  earlier  case  [d,_  ■ 10  X,  X ■ lOOOX , S/N  * 10].  Eq.  (67) 

u n 


new  yields 


D(0)  - 


sin  0 


degrees 


Comparing  with  Eq.  (54)  we  find  that  the  3 subarvay  configuration  Is  poorer 

-1/2 

by  a factor  which  Is  almost  precisely  equal  to  (1  - p . In  Eq.  (54)  we 

assumed  prior  knowledge  to  the  effect  that  the  source  was  remote  compared  with 
the  Interarray  spacing.  Eq.  (71)  makes  no  such  assumption.  Bearing  and  range 
must  therefore  be  estimated  simultaneously  with  a resultant  degradation  In  each 
estimate. 


Perhaps  the  most  Interesting  comparison  Is  that  between  the  3 subarray  system 

using  envelope  Information  only  and  a single  subarray  system  which,  because  of 

the  close  spacing  of  sensors,  can  use  the  full  narrowband  signal  without 

encountering  serious  ambiguity  problems.  The  bearing  error  of  the  single 

subarray  system  Is  given  by 

D 2/ex  . ^ -J_  1 S/N  , _ X .2 
°1  " 2ir  tS;:  sin  6^ 
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where  S/N  is  sclll  the  poHt-beamformlng  signal  to  noise  ratio.  Dividing 
(72)  by  (67)  ona  obtains 


_K2  1 + S/N 

xj  1+3  S/N 
H 


» - 


(73) 


Depending  on  various  parameter  values  this  quantity  can  be  either  larger  or 
smaller  than  unity.  Suppose,  In  particular,  that  d^^  " 100  L,  ■ lOOX  (1% 
bandwidth)  and  S/M  >>  1.  Then 

2 “ 

u/(0) 

■A-  v-  0.15  (74) 

0*^(0) 


For  this  particular  combination  of  parameters  the  single  array  is  therefore 
preferable.  In  this  Instance  the  factor  ( 1 - P^q)  plays  a key  role  In  degrading 
the  performance  of  the  multiple  array  system.  For  the  single  array,  with  Its 
relatively  short  length  L,  It  Is  quite  realistic  to  assume  that  the  source  Is  In 
the  far  field.  For  the  combination  of  subarrays  a similar  asstimptlon  might 
be  quite  unrealistic  and  one  Is  therefore  forced  Into  a joint  bearing  and  range 
estimation  problem. 

The  corresponding  computations  for  range  estimation  are  quite  straightforward. 
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required  observation  time  therefore  varies  as  (Au)  and  can  easily  become 
large  If  the  signal  bandwidth  Is  small. 

For  the  single  array  working  with  a narrowband  signal  the  mean  square  range 


error  Is 


_ 2.A.  45  1 . r X .2  1+  S/N  1 

D,  (r)  - TT  rr-  (-5 5-)  ^ — (77) 

1 8n  TAo.  ^2  ^^^2g  ^^^^2  1 - 

-1 

Eq.  (77)  contains  the  degradation  factor  (1  - because  range  measurement 

with  the  single  subarray  Is  possible  only  If  the  signal  wavefront  has  significant 
curvature  over  the  array.  When  the  source  Is  actually  remote,  so  that 
close  to  the  asymptotic  value  given  In  Eq.  (70),  then  p^^  will  certainly  be  close 
to  the  same  asymptotic  value.  It  follows  from  Ec[s.  (77)  and  (68)  that 
2 '' 

D^(r)  1 + 3 S/N  L ^ ^ 

Suppose  d^^  ■ 100  L and  consider  the  very  small  bandwidth  X^  * 1000  X.  Then  for 
S/N  » 1 


2 * 
D/(r) 


2 ■* 
D^(r) 


=>  200  . 


Thus  the  multiple  array  system  Is  clearly  preferable,  as  one  might  expect.  A 
similar  conclusion  should  hold  for  most  cases  In  which  the  subarray  spacing 
at  all  large  cosq>ared  with  the  dimensions  of  the  Individual  subarrays. 

c)  Search  area  considerations. 

The  ras  bearing  and  range  errors  Just  calculated  may  be  somewhat  misleading 
as  a measure  of  the  accuracy  with  which  the  source  can  be  located.  Fig.  12 
Illustrates  the  relevant  considerations.  We  have  previously  obtained  the  ras 


W.  J.  Bangs,  p.  104. 
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bearing  and  range  errors  Ar  and  A0  respectively.  We  therefore  know  that, 
with  reasonably  high  probability,  the  source  Is  located  In  the  curvilinear 
rectangle  of  approximate  area  4rA0Ar  shown  In  the  diagram.  However,  If  the 
range  and  bearing  errors  are  correlated,  the  high  probability  region  may  actually 
look  more  like  the  dotted  ellipse  and  Its  area  might  be  significantly  smaller. 


O 

s u.  b - 
oirro»y  3 


Fig.  12 

To  gain  soma  Insight  Into  this  question  we  tiss  the  following  generalisation 

of  the  Crasiir-Rao  bound  to  multidimensional  estimationfx 

Let  £ be  the  vector  of  estimation  errors  and  its  covariance  matrix.  Then 

for  spy  unbiased  estimate  the  concentration  ellipse 

^ H.  L.  Van  Trees,  Detection,  Estimation,  and  Modulation  Theory,  Part  1,  p.  81. 
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e A e - C 

— e — 


lies  outside  or  on  the  ellipse 


T T ^2 

e J e ■ C 


(80) 


(81) 


where  J Is  the  Fisher  Information  matrix.  As  In  the  one  dimensional  case, 
this  bound  can  be  approached  arbitrarily  closely  for  very  long  observation 
times.  We  therefore  calculate  the  area  of  the  ellipse  given  by  Eq.  (81), 
choosing  the  constant  C somewhat  arbitrarily  as  unity. ^ 

For  the  three  subarray  configuration  the  elements  of  the  Fisher  Information 
matrix  are 


J - . (a^)2  Trd^t  +). 

ij  3Tr  1 + 3 S/N 


(82) 


The  Indices  1 and  J can  assume  the  values  r and  0.  It  is  obvious  by  Inspection 


that  our  previous  results,  Eqs.  (55)  and  (56),  are  siiq>ly  the  elements  [J 

and  [J  ^1  of  the  Inverse  Fisher  Information  matrix. 

‘ rr 


ee 


The  uncertainty  ellipse  described  by  Eq.  (81)  with  C ■ 1 has  semi-4iajor  axis 
and  semi-minor  axis  X^,  given  by  the  eigenvalues  of  the  covariance  matrix 
Cov(rA9,  Ar).  Its  area  Is  therefore 


A ■ 7r/Xj^X2 

Since  the  true  range  r Is  simply  a constant,  Cov(rAG,  Ar)  Is  very  slnq>ly 
related  to  the  Fisher  Information  matrix 


(83) 


Figure  13 


V*  thereby  defining  the  region  of  high  probability  as  that  in  which  the 
exponent  of  the  (Gaussian)  error  distribution  la  no  smaller  than  -1. 
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(84) 


Eq.  (84)  defines  as  the  variance  of  (rAQ)  and  ^^2  ^ variance  of  Ar. 


The  eigenvalues  of  Cov(rA@,  Ar)  are  given  by  the  solutions  of  the  equation 

*^11^22  " ^12^  " 

It  follows  that 


^1^2  “ ^11^22  " *^12 


Det[Cov(rA0 ,Ar) ] 


Det  J * 


(86) 


A few  lines  of  algebra  now  yield  the  area  of  the  uncertainty  ellipse 


3it  (1  + 3 S/N) 


A - 


Au  T(S/N)‘^  (Ao))'^  jC~',  TZr~i  ^ /I  - p „ 

»Tr(T^  )Tr(T  ) ^rO 

The  area  of  the  rectangular  uncertainty  region  In  Fig.  12  Is  given  by 

: 1 


A " 4rA0Ar 
r 


12ti(1  + 3 S/N) 


1 - p. 


(87) 

Uance  the  ratio  of  the  two  uncertainty  areas  Is 
A 


_r 

A 


4 

Tt 


(88) 


r0 


The  factor  4/ir  la  close  to  unity  and  Is,  In  sny  case,  somewhat  arbitrary, 
since  It  depends  oii  the  choice  of  the  constant  C In  Eq,  (81).  The  significant 
conclusion  Is  chat  proper  exploitation  of  tha  dependence  between  bearing  and 
range  estimation  can  reduce  the  effective  search  area  by  a factor  proportional 
t.  (1  - 


i i 
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8.  Upper  bound  on  estimation  error.  An  envelope  tracker. 

In  the  previous  section  we  set  lower  bounds  on  the  estimation  errors 
attainable  with  any  envelope  tracker.  To  do  this  we  assumed  that  the 
pre-envelope  functions  were  available  as  separate  Inputs  to  an  optimal 
estimator.  In  actual  fact  the  pre-envelope  functions  are  not  available 
to  the  observer  and  as  a result  there  Is  no  reason  to  expect  that  the 
calculated  error  levels  are  really  at.:alnable.  In  the  present  section 
we  compute  the  estimation  errors  of  a particular  envelope  tracker.  No 
claim  Is  made  concerning  the  optimality  of  this  Instrumentation,  but  since 
It  Is  clearly  realizable,  the  calculated  errors  will  set  upper  limits  on  the 
minimum  estimation  errors  attainable  with  any  realizable  Instrumentation. 

To  keep  algebraic  complexity  within  reasonable  bounds  we  consider  a 
configuration  using  only  two  sub arrays  and  deal  with  the  estimation  of 
differential  delay.  Even  In  a multiple-array  system  all  of  the  relevant 
Information  Is  contained  In  the  differential  delay  between  the  signals 
received  by  various  subarray  p£drs.  As  far  as  the  comparison  between  the 
upper  and  lower  bounds  Is  concerned,  one  therefore  feels  that  the  solution 
for  two  subarrays  should  be  quite  representative  of  the  more  general  problem. 

The  system  to  be  analyzed  Is  shown  In  Fig.  14.  It  Is,  In  effect,  an  ordinary 
two  element  split  beam  tracker  working  with  the  square  of  the  envelope.  The 

A 

desired  estimate  ol  differential  delay  Is  the  delay  T required  In  the  lower 
channel  to  cause  the  output  z to  be  zero.  For  Gaussian  signals  (and  large  TW 
products)  this  Instnimentation  Is  known  to  be  optimal.  Its  operating  principle 
la.  In  essence,  to  ccose  correlate  the  two  channels  and  determine  the  delay  for 
which  the  cross-correlation  Is  a maximum.  Since  the  correlation  properties 
of  a random  process  are  not  changed  drastically  by  a square  law  operation,  one 
has  some  hope  that  this  procedure  may  still  be  near-optimal  for  the  non-Gausslan 
case  at  hand.  The  use  of  the  squared  envelope,  rather  than  the  envelope  Itself, 


as  the  effective  system  Input  It  simply  a matter  of  mathematical  convenience. 

2 

R can  certainly  be  constructed  from  the  data  as  readily  as  R and  the  compu- 
tation therefore  yields  a legitimate  upper  bound  on  the  minimum  attainable 
error. 


Fig.  14 

If  the  smoothing  time  T of  the  Integrator  Is  significantly 
the  envelope  correlation  time,  z(t)  Is  approximately  Gaussian. 

Is  then  unbiased  and  Its  mean  square  error  Is  approximated  by^ 

D^(z)U 


larger  than 
The  estimator 


(89) 


T Is  the  true  value  of  differential  delay  and  z Is  the  expected  value  of  z. 


W.  J.  Bangs,  S4.3.  Eq.  (89)  assumes  that  gross  ambiguities  have  been  resolved. 
When  only  envelope  Information  la  used,  this  should  not  present  a problem. 


The  denominator  of  Eq.  (89)  is  easily  evaluated 


z - y - E{~[x  ^(t)  + X ^(t)][x  ^(t  - t)  + X ^(t  - t))} 

dt  82 


(90) 


Hence 


--AE{[x  (t)x  (t)+x  (t)x  (t)][x  (t-T)x  (t-T)+x  (t-t)x  (t-i)]} 
^1  ^1  ^1  ^1  ^2  ^2  ^2  ^2 

(91) 


x(t)  stands  for  the  time  derivative  of  x(t). 

All  of  the  pre-envelope  functions  In  Eq.  (91)  are  Gaussian  random  processes. 
The  required  fourth  order  moments  are  therefore  readily  calculated.  The  result 
Is  particularly  simple  because  a zero  mean  Gaussian  process  and  Its  derivative 
are  uncorrelated,  l.e. 


E{x(t)x(t)}  - 0 


(92) 


For  narrowband  processes  with  spectra  symmetrical  about  the  center  frequency,  we 
concluded  In  section  7 that  the  pre-envelopes  x^(t)  and  Xg(t  - t)  are  uncorrelated, 
regardless  of  the  delay  T.  Finally,  If  the  noise  Is  spatially  Incoherent,  the 

A 

noise  components  contribute  nothing  to  the  average.  The  delay  t • t aligns 
the  signal  components  In  the  two  channels  and  one  Is  left  with 


dz 

dt 


- -8  x^  ^(t)  • i ^(t) 

Cl 


(93) 
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If  the  signal  preoenvelope  has  autocorrelation  R (t),  this  reduces  to 

s 


dz 

di 


- 8 R (0)  R "(0) 

a s 


(94) 


T«T 


where  R "(0)  is  the  second  derivative  of  the  signal  autocorrelation  evaluated 


at  the  origin 
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The  numerator  of  Eq.  (89)  Is  a good  deal  more  tedious  to  compute.  The 
procedure  Is  straightforward  In  principle.  The  variance  of  z Is  simply  the 
Integral  of  Its  spectrum  G^(u)). 

00 

2 ( 

D (z)  - G (w)  do)  (95) 

Jo  * 

The  spectrum  of  z Is  readily  obtained  from  the  spectrum  of  y [see  Fig.  14]. 

G (w)  - G (w)  |H(a))|^  (96) 

z y 

where  U(u)  Is  the  transfer  function  of  the  Integrator.  G^(ii)),  In  turn.  Is 
computed  from  the  Fourier  transform  of  the  autocorrelation  R (a). 

00  y 

Gy(u)  - 7 I Ry(“)e"^““  da  . (97) 

—00 

The  computational  problem  therefore  revolves  about  the  determination  of 
Ry(o)  - E{y(t)  y(t  + a) } 

- 4E{[x  (t)x  (t)  + X (t)x  (t)][x  ^(t  - t)  + X ^(t  - t)] 

1 *^1  ®1  ®1  ^2  ®2 

• [x  (t+a)x  (t+a)  + X (t+a)x  (t+a)](x  ^(t-T+o)  + x ^(t-x+a)]}  (98) 

®1  *^1  ®1  ®1  *^2  ®2 

Multiplying  out  the  four  binomials  one  obtains  a number  of  terms  each  of  which  *• 

Is  an  8^^  order  moment  of  Jointly  Gaussian  random  variables.  Each  such  8**^ 

order  moment  consists  of  all  possible  combinations  of  second  order  moments 

V taken  four  at  a time,  of  which  there  Is  a very  large  number.  Many  of  the 

terms  are  Identical  and  the  computational  problem  Is  almost  entirely  one  of 

keeping  an  accurate  count  of  terms  of  various  types.  Fortunately  there  is 

a svetemetic  scheme  for  doing  thls^.  The  algebraic  detalla  are  extremely 

^ J.  H.  Lanlng  and  R.  H.  Battin,  Random  Processes  in  Automatic  Control, 

McGraw-Hill,  1956,  p.  82-85. 
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tedious  and  without  Interest  In  this  discussion.  We  therefore  only  state  the 
final  result 

Ry(a)  - 8{[-R^(a)R'^(a)-(R’(a))^][16R^2(0)+4Rjj^(0)+16R^(0)Rjj(0) 

+ 8R^^(a)  + 8Rjj^(a)  + 8R^(o)Rjj(a)  ] 

+ [-Rjj(a)R^(o)  - RjJ(a)R^(o)][8R^(0)  + 4Rjj^(0)  + 8R^(0)Rjj(0) 

+ 8R^^a)  + 4Rj^^(o)  + 4Rjj(o)R^(o)] 

+ [-R^(a)Rjj"(c)  - R'(o)Rj^'(o)][8R^^(0)  + 4Rjj^(0)  + 16  R^(0)Rjj(0) 

+ 8R^^(a)  + 4R^^(a)  + IZR^j (o) R^ (a)  ] 

+ [-Rj,(a)R;j(a)  - [4R^2(0)  + 4Rjj^(0)  + 8Rjj(0)R^(0) 

+ 4R^^(o)  + 4Rjj^(o)  + 8R^(a)Rjj(a)  ] 

+ 16(R'(a))^[R^^(o)  + R^(a)Rjj(a)]  - 8R^(a)R^(o)R^(o)Rjj(o) 

- 8R^(a)R;|(a)R^^0)  - 4Rjj(o)R;J(o)R^^(0)  }.  (99) 

i^(a)  la  the  autocorrelation  of  the  noise  pre-envelope. 

To  proceed  further  we  must  assume  a particular  form  of  the  signal  and  noise 
pre-envelope  autocorrelations.  For  computational  convenience  we  choose 


2 2 

0 a 

R (o)  - Pe“  2 

8 

(100) 

1.  2 2 
ko  a 

Rj^(o)  - Ne  2 

(101) 

P la  the  total  signal  power  and  N the  total  noise  power,  k is  a numerical 
constant  which  governs  the  relation  between  the  signal  and  noise  bandwldths. 

In  practice  the  noise  bandwidth  would  be  determined  primarily  by  the  bandwidth 


of  the  processing  system  (since  the  received  noise  Is  likely  to  have  a spectrum 
flat  over  a band  much  larger  than  that  occupied  by  the  signal) . The  processing 
band  could  at  best  be  matched  to  the  signal,  leading  to  a value  k ~ 2 (because 
the  processing  filter  also  affects  the  signal).  If  the  center  frequency  of 
the  narrowband  signal  Is  not  known  precisely,  one  might  use  a substantially  larger 
processing  bandwidth  and  hence  a larger  value  k. 

Substituting  Eqs.  (100)  and  (101)  Into  Eq.  (99)  and  evaluating  the  Fourier 
Integral  of  Eq.  (97)  one  obtains 
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(102) 
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£q.  (102)  Is  still  very  cumbersome.  It  Is  reproduced  in  full  primarily  in 

order  to  make  the  following  observation.  The  equation  has  two  types  of  terms: 

2 2 

Those  which  contain  a factor  (u  /a  ) and  those  which  do  not.  There  are  14 

of  the  former  and  only  6 of  the  latter,  a Is  the  bandwidth  of  the  signal. 

If  To  >>  1 (large  TW  product)  the  bandwidth  of  the  T second  Integrator  In 

Fig.  14  Is  small  conq)ared  with  o,  so  that  the  terms  containing  the  factor 
2 2 

(u  /a  ) contribute  little  to  the  output  variance  compared  with  terms  which 

2 2 

have  no  such  factor.  Before  discarding  the  (u)  /a  ) terms,  however,  It  Is 

4 2 2 

Important  to  note  that  there  are  no  terms  In  P which  do  not  contain  (u)  /c  ). 

2 2 

To  discard  all  terms  containing  (u  /a  ) would  lead  to  the  obviously  false 

result  that  the  output  variance  Is  zero  In  the  absence  of  noise.  More 

2 2 

specifically,  we  shall  find  that  the  terms  not  containing  (m  /o  ) yield 

contributions  to  the  output  variance  of  order  (To)  while  those  containing 

2 2 -2 
(u  /o  ) yield  contributions  of  order  (To)  . For  large  (To)  products  and 

moderate  signal  to  noise  ratios  the  former  are  clearly  dominant.  For  large 

signal  to  noise  ratios  and  only  moderate  (To)  products  the  latter  may  assume 

primary  Importance. 

There  remains  the  Identification  of  the  Integrator  transfer  function  H(u) 
and  the  evaluation  of  the  output  variance  from  Eq.  (95).  The  transfer 
function  of  a T second  integrator  Is 


. . 1 1 - 
H(a,)  — 


(103) 


The  quantity  |H(b>)|  required  In  Eq.  (96)  Is  therefore 
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pass  loaedlately  to  Che  final  result  for  estimation  error. 
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(105) 


Mien  To  » 1 the  second  term  In  each  power  of  N/P  would  tend  to  be  small  compared 

with  the  first.  As  anticipated,  however,  the  (N/P)*  contribution  only  has  a 
-2 

term  of  order  (To)  . For  very  large  signal  to  noise  ratios  the  magnitude 
of  Che  To  product  therefore  becomes  critical. 

To  interpret  Che  k dependence  of  Eq.  (105)  properly,  it  is  necessairy  to 
keep  in  mind  that  N is  the  total  noise  power.  If  Che  underlying  noise  is  white 
and  its  effective  bandwidth  is  determined  by  the  processor  bandwidth,  Chen  N 
is  proportional  to  i^T.  In  chat  case  the  M/P  term  is  asymptotically  independent 
of  k,  all  other  noise  dependent  terms  grow  with  k. 
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9.  Comparison  of  upper  and  lower  bounds 

The  primary  purpose  of  analyzing  the  split  beam  envelope  tracker  for  the 
simple  configuration  of  two  subarrays  was  to  check  whether  It  could  estimate 
relative  delay  with  anything  like  optimal  efficiency.  We  must  therefore 
compare  Eq.  (105)  with  the  comparable  result  In  section  7a.  Unfortunately 
the  results  In  the  two  sections  are  not  directly  comparable  as  they  stand. 

In  order  to  keep  down  computational  complexity  we  used  the  correlation 
functions  given  by  Eqs.  (100)  and  (101)  [corresponding  to  Gaussian  spectra] 

In  section  8.  In  order  to  avoid  other  computational  difficulties  in  section 
7a  we  worked  with  signal  and  noise  spectra  flat  (or  at  least  of  identical 
shape)  over  the  signal  band.  The  result  In  section  7a  [Eq.  (SI)]  therefore 
Involves  the  signal  bandwidth  parameter  Au  which  is  not  simply  related  to 
the  equivalent  parameter  a In  Eq.  (105).  To  make  the  upper  and  lower  bounds 
comparable  we  return  to  section  7a,  use  the  spectra  Introduced  In  section  8 
and  circumvent  the  computational  difficulty  by  finding  a lower  bound  on  the 
resulting  Integral.  We  are  therefore  obtaining  a lower  bound  on  the  lower  bound 
which,  together  with  Eq.  (105),  still  successfully  brackets  the  attainable 
minimum  mean  square  error. 

The  spectra  corresponding  to  Eqs.  (100)  and  (101)  are 
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(106) 


N(u)"/^  — ^ e 2ka^ 

/ka 


(107) 


Substituting  these  spectra  Into  Eq.  (50)  we  obtain 
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The  upper  bound  le  generated  by  omitting  the  first  denominator  term  for  0 ^ u £ 
and  the  second  denominator  term  for  u > u^.  The  frequency  will  be  chosen 
so  that  It  minimizes  the  right  side  of  Eq.  (108) , thus  setting  the  tightest 
possible  bound.  The  last  version  of  Eq.  (108)  Is  readily  Integrated 
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The  error  function  erf(x)  la  defined  by 


erf(x) 


X - 

i I - • 


(110) 


To  find  the  best  value  of  we  differentiate  the  last  version  of  Eq.  (108) 


and  set  the  derivative  equal  to  zero 
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A few  steps  of  algebra  yield 

r\  ? 

U)i  - 0 / £n  (4k  -j)  . (112) 

Not  surprisingly,  Eq.  (112)  simply  gives  the  value  of  u at  which  the  second 

denominator  term  in  the  exact  version  of  Eq.  (108)  becomes  equal  to  unity. 

2 2 

If  4k  P /M  < 1,  Eq.  (112)  is  meaningless  and  one  uses  ■ 0.  Again,  this  is 
not  surprising,  because  the  signal  to  noise  ratio  (after  combination  of  the 
subarray  outputs)  is  now  everywhere  smaller  than  unity. 

Substituting  Eq.  (112)  into  Eq.  (109)  one  obtains  after  some  algebraic 
simplification 


tD^(T)r^6-^  a^(T-~)^^^  {A:|  erf[£n(2/k  |)]^^^ 

+ — k ^ [1  - erf  (^n  4k  ^)^^^] 

/l  N N 


- ~ / ^n(2/1^  h } 


(113) 


In  practice  the  noise  bandwidth  is  almost  certain  to  be  larger  than  the 
signal  bandwidth  (k  > 1) . In  that  case  the  error  functions  very  rapidly  approach 
unity  as  tbs  signal  to  noise  ratio  F/N  becomes  larger  than  one.  Eq.  (113) 
is  than  dominated  by  the  flrat  term  and  one  can  write  to  an  excellent  approximation 

a 

' ' - 2 To  P .2 

o k 

For  vary  large  Ib  and  high  signal  to  noise  ratio  the  split  beam  tracker 
error  [Eq.  (105)]  can  be  approximated  by 

^ This  expression  properly  goes  to  saro  for  k ■ 1,  because  the  signal  to  noise 
ratio  is  Chan  constsnc  at  all  fraquancias. 


(114)- 
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(115) 


Eqs.  (114)  and  (115)  differ  only  in  the  k dependent  term.  In  this  limit  of 
large  To  and  high  signal  to  noise  ratio  one  can  therefore  write  the  ratio  of 
the  upper  bound  to  the  lower  bound  as 


2 

Upper  bound  k , k , 3k  1 , zntx 

Lower  bound  " . i)3/2  + i)3/2  + 3)3/2  * 

Over  the  practically  significant  range  of  k (k  ^ 2)  Eq,  (116)  Is  surprisingly 
Insensitive  to  k.  It  assumes  a value  of  approximately  4 at  k « 2,  declines 
to  a minimum  of  about  3.25  near  k ■ 5 and  then  rises  to  an  asymptotic  value 
of  4 as  k -»■  <».  Thus  the  lower  bound  falls  5 to  6 db  below  the  upper  bound, 
essentially  Independent  of  k. 

To  Interpret  this  result  properly  one  must  keep  In  mind  that  the  lower 
bound  of  Eq.  (50)  Is  almost  certainly  unattainable  because  It  assumes  the 
separate  availability  of  the  pre-envelope  functions.  Furthermore,  Eq.  (113) 
only  computes  a lower  bound  on  this  lower  bound.  It  would  therefore  not  be 
surprising  to  find  Eq.  (113)  low  by  as  much  as  3 db.  If  this  Is  the  case, 
there  is  at  least  a range  of  To  and  P/N  for  which  the  simple  split  beam  envelope 
tracker  comes  within  3 db  of  the  absolute  optimum.  One  must  further  keep  In 
mind  that  no  effort  has  been  made  to  compensate  the  tracker  In  Fig.  14  for 
signal  and  noise  spectral  properties.  A filter  of  the  Eckert  type  in  each 
channel  would  very  probably  Improve  Its  performance  relative  to  the  lower 
bound,  quite  possibly  leaving  a differential  too  small  to  be  of  practical 
concern. 

The  discrepancy  between  the  upper  and  lower  bound  becomes  much  more 
pronounced  when  either  of  the  assumptions  To  »>  1 and  F/H  » 1 falls. 

Consider  first  the  case  of  small  signal  to  noise  ratio. 

^ ^ JT-  ■ if  ■ 

« - — - - - — - - ^ ' ■ . 
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^or  P/N  <<  1,  £q.  (103)  Is  domlnsted  by  the  terms  In  (N/P)  and  (N/P) 
which  are  not  present  In  the  lower  bound  at  all.  The  square  law  envelope 
detector  and  multiplier  together  generate  noise  intermodulation  products 
of  order  as  high  as  4.  When  P/N  « 1 these  contribute  most  of  the  output 
fluctuation  of  the  split  beam  envelope  tracker.  For  low  signal  to  noise 
ratios  one  therefore  has  strong  reasons  to  suspect  that  the  simple  instrumen- 
tation of  Fig.  14  falls  well  short  of  the  optimum.  In  practice,  of  course,  no 

m 

tracker  will  perform  well  unless  P/N  (the  post-beamf ormlng  signal  to  noise  ratio 
at  each  subarray)  is  reasonably  large.  This  particular  shortcoming  may  therefore 
be  a somewhat  academic  matter. 

A performance  limitation  of  at  least  potential  practical  Impact  results 
from  the  tracker's  behavior  at  moderate  To  products.  If  the  signal  bandwidth 
is  very  small  it  may  not  be  possible  to  choose  T so  large  that  To  exceeds 
unity  by  many  orders  of  magnitude.  In  such  a situation  the  first  (noise 
independent)  term  of  Eq.  (105)  may  make  the  greatest  contribution  to  the 
estimation  error.  Once  this  happens,  further  Increases  in  signal  to  noise 
ratio  will  no  longer  have  much  effect  on  the  error,  which  is  now  determined 
primarily  by  signal  fluctuations.  A plot  of  mean  square  error  versus  P/N 
will  therefore  approach  a minimum  given  by 


2 * 
D^(t) 


min 


5 _1 


(To)^ 


(117) 


The  lower  bound  of  Eq.  (113),  by  contrast,  approaches  zero  as  P/N  -*■  <». 

The  signal  to  noise  ratio  at  which  Eq.  (105)  approaches  Eq.  (117)  depends 
on  Che  relative  magnitude  of  the  coefficients  of  (N/P)^  end  (N/P)^.  For  k 
in  the  range  2 ^ k ^ 10  one  finds 
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To. 


(119) 


If  To  is  only  of  the  order  of  10,  signal  to  noise  ratios  satisfying  Eq.  (119) 


might  well  be  encountered  in  practice. 

Actual  plots  of  mean  square  error  as  a function  of  signal  to  noise 
ratio  are  given  in  Figs.  15-17  for  several  combinations  of  parameter  values. 

Fig.  15  uses  To  ■ 10,  k - 2;  in  Fig.  16  To  - 100,  k - 2.  As  anticipated, 
the  upper  bound  in  Fig.  15  approaches  its  minimum  near  P/N  - 10  db,  whereas 
in  Fig.  16  it  Only  levels  off  near  P/N  ■ 20  db.  Fig.  17  uses  To  ■ 100,  k ■ 10. 

Comparison  with  Fig.  16  indicates  that  the  effect  of  the  k variation  is  j 

I 

relatively  minor.  In  each  case  the  upper  bound  approaches  to  within  about 
6 db  of  the  lower  bound  befor  it  levels  off  at  very  high  signal  to  noise 
ratios.  At  signal  to  noise  ratios  near  0 db  one  sees  the  effect  of  the 
(N/F)^  terms  in  Eqs.  (105)  and  (113).  The  region  P/N  < 0 db  was  omitted 


P/N  (db) 
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10.  An  Improved  Envelope  Tracker 
In  Che  previous  section  we  have  seen  that  Che  simple  split  beam 
envelope  tracker  falls  to  approach  Che  theoretical  optimum  at  every  high 
signal  to  noise  ratios  unless  the  Clme-bandwldCh  product  Is  extremely  large. 
Since  the  condilnatlon  of  high  signal  Co  noise  ratios  and  only  moderately 
large  tlme-bandwldth  products  Is  not  at  all  unlikely  In  practice , we  look 
for  an  Improved  Instrumentation. 

On  a purely  formal  level,  we  note  that  the  difficulty  can  be  traced 

to  the  first  (noise  Independent)  term  In  Eq.  (105).  This  component  of  Che 

mean  square  error  Is  derived  from  the  first  term  of  Eq.  (102)  by  use  of 

Eqs.  (96)  and  (95).  In  our  case  the  low  pass  filter  Is  a simple  T second 

2 

Integrator  with  the  transfer  function  H(u)  given  by  Eq.  (103).  |H(b))| 

decays  only  with  the  second  power  of  frequency.  As  a result  terms  of  Che 
2 2 

form  0)  exp(-au  ) In  Eq.  (102)  make  substantial  contributions  Co  the  mean 
square  error.  In  particular,  the  noise  Independent  term  Is  of  this  form  and 
can  therefore  not  be  Ignored. 

This  line  of  reasoning  Immediately  suggests  a modification  of  the  tracker. 
One  need  only  replace  Che  T second  Integrator  with  a low  pass  filter  whose 
transfer  function  decays  more  rapidly  with  frequency.  Computations  were  carried 
out  for  Che  low  pass  filter  2 

H(u))  - e ^ 0 . (120) 

The  frequency  ratio  o/Oq  now  takes  Che  place  of  the  tlme-bandwldth  product 
To.  One  finds  that  the  curves  of  mean  square  error  versus  signal  to  noise 
ratio  behave  qualitatively  like  the  upper  band  curves  of  Figs.  15-17,  but 

that  the  final  level  Is  reached  at  a signal  Co  noise  ratio  of  the  order 

2 -2 
(o/Oq)  • The  magnitude  of  the  final  mean  square  error  Is  of  order  (o/Oq) 

rachar  than  (To)  A definite  Improvement  has  therefore  been  made. 
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We  shall  not  pursue  the  subject  of  modifying  the  low  pass  filter  because 
It  Is  possible  to  take  a more  fundamental  approach.  The  troublesome  term 
In  Eq.  (105)  describes  the  fluctuations  produced  by  Intermodulatlng  signal 

A 

components.  When  the  delay  t in  Fig.  14  Is  adjusted  properly  the  signal 
components  from  the  two  subarrays  are  Identical  and  It  should  be  possible  to 
cancel  their  contribution  to  the  output  fluctuation.  This  Idea  Is  exploited 
In  the  Instrumentation  shown  In  Fig.  18.  We  shall  refer  to  It  as  a balanced 
split  beam  envelope  tracker.^  It  consists  of  2 trackers  using  the  same  data 
whose  outputs  are  subtracted.  In  one  of  the  units  the  output  of  subarray  1 
Is  differentiated.  In  the  other  that  of  subarray  2.  The  Input  Into  the 
Integrator  Is  therefore  given  by 


y(t)  - R^^(t)]  R2^t  - X)  - Ri2(t)[^  R2^(t  - x)  ] 


f_jd  ^1 

^dt  _ 2 


R2"(t  - T) 


] R2  (t  - x). 


(121) 


When  there  la  no  noise  and  x Is  adjusted  to  align  the  signal  components, 
R^^(t)/R2^(t  - x)  - 1 and  the  right  side  of  Eq.  (121)  Is  Identically  equal 
to  zero.  Thus  the  pure  signal  contributions  to  the  output  fluctuation  are 
eliminated  at  the  operating  null.  There  will,  of  course,  be  signal  components 

A 

In  the  error  trtien  x Is  not  perfectly  adjusted. 

The  analysis  of  Fig.  18  proceeds  exactly  as  that  of  Fig.  14.  For 
Chat  reason  we  only  state  the  final  result 


The  ides  of  balancing  Che  tracker  configuration  and  thereby  removing  signal 
contributions  to  the  output  fluctuation  is  equally  applicable  when  the  data 
vector  consists  of  the  complata  sensor  outputs  rather  than  their  envelopes.  The 
usual  (unbalance<9  tracker  is  optimal  only  under  the  asaumption  TV  »’  •*. 
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(122) 


ComparlsoQ  with  Eq.  (105)  reveals  the  following  features: 

1)  Eq.  (122)  does  not  contain  a noise  Independent  term.  This  Is,  of 
course,  the  feature  which  motivated  construction  of  the  balanced  tracker. 

2)  The  coefficient  of  each  power  of  (N/F)  In  Eq.  (122)  Is  smaller  than 
the  corresponding  coefficient  of  Eq.  (105).  This  Is  true  for  all  values 
of  k and  lb . Hence  the  balanced  tracker  Is  better  than  the  conventional 
tracker  under  all  operating  condition. 

3)  For  reasonably  large  To  the  coefficient  of  the  N/F  term  In  Eq.  (122) 

Is  only  slightly  smaller  than  that  of  the  corresponding  term  In  Eq.  (105). 
In  this  Important  region,  therefore,  the  balanced  tracker  Is  only 
marginally  better  than  the  conventional  version. 

A typical  plot  of  mean  square  error  versus  signal  to  noise  ratio  for  the 
balanced  tracker  Is  given  In  Fig.  19.  Also  shown  for  comparison  Is  the  lower 
bound  and  the  performance  curve  for  the  coavantloual  (rocker.  The  parameter 

values  are  To  • 10,  k ■ 2.  The  curve  for  the  balanced  tracker  differs  from 
the  lower  bound  by  little  more  than  a vertical  offset  of  about  6 db.  Fotentlal 
Improvamenta  are  therefore  limited  to  an  absolute  maximum  of  6 db.  In  actual 
fact  tha  room  for  Improvemant  la  probably  smallar  bacause  tha  lower  bound  Is 
not  realisable  [aee  discussion  In  eectlon  9].  The  balanced  split  beam  tracker, 
possibly  Improved  by  the  Introduction  of  Eckart  flltars,  may  therefore  be 
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a very  reasonable  approximation  to  the  optimum. 

Further  work  on  the  balanced  tracker  la  planned  for  the  future.  Aside 
from  the  effect  of  Eckart  filters  and  similar  modifications  one  needs  to 

A 

know  the  output  fluctuation  when  t Is  not  perfectly  adjusted.  Under  high 
signal  to  noise  ratio  conditions  one  would  expect  the  output  fluctuations  to 

A 

exhibit  a sharp  dip  near  t ■ t.  This  effect  may  actually  provide  a useful  tool 

A 

for  Identifying  the  proper  value  of  t. 


) 


